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Irreversible thermodynamics of single-molecule experiments subject to external constraining forces 
of a mechanical nature is presented. Extending Onsager’s formalism to the non-linear case of systems 
under non-equilibrium external constraints, we are able to calculate the entropy production and the 
general non-linear kinetic equations for the variables involved. In particular, we analyze the case 
of RNA stretching protocols obtaining critical oscillations between different configurational states 
when forced by external means to remain in the unstable region of its free-energy landscape, as 
observed in experiments. We also calculate the entropy produced during these hopping events, and 
show how resonant phenomena in stretching experiments of single RNA macromolecules may arise. 
We also calculate the hopping rates using Kramer’s approach obtaining a good comparison with 
experiments. 
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I. INTRODUCTION 

Many mesoscopic systems of great interest ranging 
from biology, such as proteins, RNA or DNA molecules 
[IHIO], to nano-engeneering, like rotaxanes and catenanes 
[IIHISl, have similar dynamic behaviors due to their in¬ 
trinsic properties and structure, as well as to their strong 
coupling with the surroundings. In particular, when they 
are subject to stressing conditions, like stretching forces 
or radiative interaction, the dynamics of these systems 
may present critical oscillations and stochastic resonance 
phenomena. Although well understood from the point 
of view of stochastic processes theory, the connection be¬ 
tween these dynamical behaviors and the thermodynamic 
properties of these small systems and its coupling with 
changing surroundings is poorly understood m- 

A better understanding of these peculiarities is of great 
importance in molecular biology, since it opens the possi¬ 
bility of correlating relevant physical-chemistry laws with 
simple models, thus allowing a transparent interpretation 
of experimental data and an agile prediction of different 
effects. Here, we establish this correlation through an ir¬ 
reversible thermodynamic description of biological small 
systems subject to external forces of a mechanical nature. 
This formalism is a powerful generalization of a previous 
analysis performed in Ref. |15j involving thermal con¬ 
straints and permits one to calculate non-linear kinetic 
equations for the variables comprised and the entropy 
produced during the hopping events in both the critical 
oscillations and stochastic resonant phenomena. In this 
way, we can better distinguish between reversible and ir¬ 
reversible behaviors of small systems and their coupling 
with the surroundings, through thermal and non-thermal 
fluctuations. 

In order to highlight the ideas and simplify the presen¬ 
tation, we apply our formalism to analyze single RNA 
macromolecule stretching experiments which provide in¬ 
formation on its energy landscape [SHU US] . Specifically, 
the existence of a loop in the force-extension curve (FEC) 


implies that there is a range of tensions for which two 
stable configurations of the system are possible. We 
will show that the transition between these two possible 
(folded/unfolded) configurations of RNA is concomitant 
with an irreversible generation of entropy by the molecule 

[IZlIIHj- 

Typically, the theoretical modeling of RNA macro¬ 
molecule stretching experiments is performed by consid¬ 
ering the composite system: RNA-macromolecule plus 
heat-bath, as an isolated system in thermal equilibrium. 
However, we propose here a novel approach to the prob¬ 
lem of single macromolecule stretching experiments that 
differs from previous theoretical analysis because it is 
based on Onsager’s non-equilibrium thermodynamics, 
specifically applied to the description of the time evo¬ 
lution of the state of the single macromolecule. This 
huge difference with respect to the mentioned previous 
descriptions, allows one to understand in a simpler and 
more consistent way the effect of fluctuations on both, 
the single macromolecule dynamics realizations as well 
as over the time averaged (or collective) properties that 
may be extracted from these experiments. Using the On¬ 
sager’s theory to obtain the evolution equation for the 
extension of a single RNA molecule is a very general and 
simple form that does not needs to consider confusing or 
even inappropriate assumptions like ergodicity and equi¬ 
librium that are commonly used in the literature on this 
problem. In general terms, the idea is to calculate the 
entropy produced by the single RNA molecule during its 
transition between the folded and unfolded states, and 
to relate this with the corresponding free energy change 
between these configurations. This procedure is very rich 
since shows that the folding/unfolding process may not 
be viewed as an equilibrium situation and permits to 
deduce the stochastic non-linear kinetic equations that 
govern the hopping process for different conditions of the 
thermal bath and of the experimental devices. 

It is convenient to mention here that some previous 
works have shown that entropy production may be as- 
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sociated to the process of transition between folded and 
unfolded states in single RNA macromolecule stretching 
experiments [niiiH]. However, they have been focused to 
show that non-equilibrium statistical arguments lead to 
the well established master equation formalism, which is 
the usual mathematical tool for the theoretical descrip¬ 
tion of these systems HU. The important precedent here 
is that it was shown that the master equation formal¬ 
ism is compatible with the fact that entropy is produced 
during the hopping dynamics of the single RNA molecule 
stretching experiments. 

The use of a master equation comes from the fact 
that it allows to determine the transition rates between 
the folded-unfolded states of the single RNA molecule, 
and suggests a chemical-like interpretation of the pro¬ 
cess. This description is sufficiently accurate from the 
quantitative point of view because it can be associated 
to average values in the thermal noise case. However, 
although non-Markovian master equations exist, using 
them to obtain practical results involving colored noises 
becomes a very difficult task that may substantially mod¬ 
ify the obtained results and even their interpretation. 
Furthermore, it may result in confusion when used under 
the perspective of chemical kinetics, as we will discuss in 
detail. The evolution equation for the extension we de¬ 
rive using our generalization of Onsager’s theory is of the 
Langevin type, and therefore can be easily generalized to 
more complex situations. 

The article is organized as follows. First, we present 
our approach to calculate the entropy production of the 
system from which we obtain and analyze the non-linear 
stochastic kinetic equation accounting for the evolution of 
the variable determining the state of the system, the elon¬ 
gation of a RNA molecule. Afterwards, we solve numeri¬ 
cally the equation showing four main results. We also cal¬ 
culate the entropy produced during a single realization of 
hopping events. Then, we discuss the ensemble descrip¬ 
tion that is widely used to reconstruct the free-energy 
landscapes of biological small-systems and discuss some 
important aspects in relation with the single molecule 
description. Finally, we present our conclusions. 


II. DISSIPATION IN SINGLE MOLECULE 
EXPERIMENTS, ENTROPY PRODUCTION AND 
KINETIC EQUATIONS 

We begin our approach to the problem by consider¬ 
ing a single RNA molecule of length (x) subject to a 
stretching force r which can be controlled as an external 
parameter. According to thermodynamics [191120) . for 
the case of constant volume and temperature, the differ¬ 
ential amount of work done by the surroundings on the 
molecule in order to change its elongation, dW = Td{x), 
satisfies the inequality 

( 1 ) 


where dF is the differential of the Helmholtz free-energy. 
As it is well known, the minimum work is performed in 
reversible processes for which the equality holds, imply¬ 
ing that dF = dWrev = Td{x). 

For unequilibrated irreversible processes, the same 
change of state comprises a larger work due to dissipa¬ 
tion and the inequality holds. Following Ref. [21], it is 
easy to show that in the irreversible case the equality can 
be restored by adding the entropy produced during the 
process, TdiS. This leads to the relation 

dW = dF + Td,S. (2) 

This equation is very general since it is valid for vari¬ 
ables determining the global state of the system (even if 
it is small) and may be used to determine the evolution 
equations of the system during the unequilibrated pro¬ 
cess. From Eq. ([^ it follows that the differential of the 
dissipated work is dWdis = TdiS. 

In order to analyze the dissipative dynamics of the 
elongation of the RNA molecule, we have to notice that 
Eq. ([^ gives the entropy produced by the system dur¬ 
ing the unequilibrated elongation process in terms of its 
free-energy change and the work done on the system 

TdiS = —dF + Td{x) . (3) 

Erom this equation the Onsager’s theory may be gener¬ 
alized by assuming the existence of linear laws between 
fluxes and thermodynamic forces. As we will show, the 
flows are given by the time derivative of the independent 
variables, d{x)/dt, whereas the forces are given as deriva¬ 
tives of the free energy with respect to the independent 
variables: dF/dx, [22]. In this sense, the non-linear char¬ 
acter of the resulting evolution equations comes from the 
fact that the free-energy is non-linear in the independent 
variables x, in contrast with original Onsager’s theory, 
in which the forces are linear functions of the fluctuat¬ 
ing variables. Finally, it is important to mention that 
in a previous work m we discussed how equations of 
the form ([^ can be derived based on statistical mechan¬ 
ics and discuss the relation of these equations with the 
internal stability of small systems. As we have shown, 
a change of convexity in the corresponding thermody¬ 
namic potential of a small system, irrespective of its na¬ 
ture, triggers irreversible processes that dissipate energy 
and, more importantly, establishes differences between 
the performed work and the corresponding free energy 
change as we discuss in the next subsections. 

A. Entropy production and stochastic kinetic 
equations 

In the general case when dF fy Td{x), a dynamics is 
established in the system. This dynamics can be ana¬ 
lyzed by computing the entropy production and using 
the second law of thermodynamics within the framework 
of Onsager’s irreversible thermodynamics jujus]. For 


dW > dF, 
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this, we may divide Eq. ([^ by dt and take the limit 
dt —>■ 0 yielding the differential equation for the entropy 
production 

dt 

Since the second law imposes that the entropy produc¬ 
tion must be non-negative, diS/dt > 0, we may assume, 
without loss of generality, the relationship 


dFi 


(x) 


d{x) 


+ T 


d{x) 

dt 


( 4 ) 


d{x) 

~dr ^ 


a(cr) 



( 5 ) 


where 7 > 0 is an Onsager coefficient having dimensions 
of mobility. It is important to mention here that an equa¬ 
tion similar to Eq. ([^ has been previously derived in Ref. 
m by using using mesoscopic non-equilibrium thermo¬ 
dynamics. Such a connection with the statistical Gibbs 
entropy is relevant for two reasons: i) It shows that dur¬ 
ing the transitions of the single molecule between config¬ 
urations entropy is produced and, ii) at the same time, 
establishes a connection between this entropy production 
and the master equation formalism mentioned in the in¬ 
troduction, see Refs. [laiiH]. In addition, it must be 
stressed that the derivation given in Ref. m introduces 
an internal degree of freedom that allows one to describe 
in detail the distribution of tension and forces along the 
system length. When this internal variable is considered, 
a thermodynamic formulation of the problem needs the 
assumption of local equilibrium in the corresponding level 
of description. The equivalence of results shows that the 
local equilibrium assumption remains valid for large gra¬ 
dients in small systems, like proteins |211 and carbon nan¬ 
otubes [25] where the non-equilibrium thermal relaxation 
as been successfully studied. 

Equation ([^ is the overdamped dynamic equation for 
the average elongation of the molecule and thus, meaning 
that at equilibrium {d{x)/dt = 0 ) the external tension r is 
compensated by the internal tension: Tint = / d{x). 

In the case of small systems in contact with a heat 
bath, like in the single RNA-molecule experiments under 
consideration, thermal fluctuations should produce fluc¬ 
tuations of the instantaneous elongation of the molecule, 
x{t), that in turn induce an unbalance of forces, dF 7 ^ 
rdx. This is the well known program of Onsager’s irre¬ 
versible thermodynamics HOlllS]. Therefore, for small 
systems one has to formulate the stochastic version of 
Eq. ([^, that is given by the Langevin-like equation 

( 6 ) 

where ^(t) is an additive random Gaussian noise having 
zero mean and obeying the fluctuation-dissipation theo¬ 
rem: {^{t)^{t')) = 2DS{t — t'), with D = the inten¬ 

sity of noise. The derivation of Eq. (|^ can also be ob¬ 
tained by applying a fluctuating hydrodynamics scheme. 

The key point here is that thermal fluctuations sponta¬ 
neously modify the state of the system which develops a 


dx 

dt 



dF^ 

dx 


+ T 


+ 



FIG. 1. (Color online) Elongation of RNA molecules subject to 
tension near the folding/unfolding states. It is also shown the corre¬ 
sponding shape of the free energy as obtained by fitting simulation 
data from [ 6 ]. a) For low external tensions applied on the RNA 
molecule the fitting of the free energy implies that the secod term 
at the right hand side of Eq. jsi is: t/// — t = —4.72pN, where 
the absolute minimum of the free energy corresponds to the folded 
configuration, b) For external tensions in the range of the unstable 
region of the force-extension curve, the free energy is bistable and 
the effective tension takes the value rj// — t = —5.24pN. c) For 
(c) 

large enough tensions, — t = —6.42pN, the absolute minimum 
of the free energy corresponds to the unfolded configuration, d) 
Effect of a periodically driving force rft) on the hopping dynamics 
of the molecule. The parameters of the oscillation are ui = 7 /IO 
and t/// — t ~ —5.32pN. 


dynamics that depends on the stochastic nature of those 
fluctuations and on the internal and external constraints. 
In the case of stretching experiments with single RNA 
molecules that we are considering, two outputs are pos¬ 
sible: If the external tension is below or above of the 

bistable region of the force-extension curve, the elonga¬ 
tion is single valued which means that these fluctuations 
should be Gaussian with respect to an average value of 
the elongation because the molecule is in equilibrium, ii) 
However, if the tension takes values in the portion of the 
force-extension curve with negative slope (corresponding 
to the convex part of the free energy) the system is intrin¬ 
sically unstable and therefore thermal fluctuations should 
trigger the critical oscillations of the elongation of the 
molecule between the folded and unfolded states. The 
probability distribution associated to these fluctuations, 
although may take a stationary shape, is bimodal and it 
cannot be taken as an equilibrium distribution compati¬ 
ble with thermodynamic equilibrium. 

According to the experimental and simulation results 
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reported in [ 3 ], for tensions around 15pN, the intrinsic 
free-energy of the RNA molecule is bistable as a function 
of the elongation x. A series expansion of around 
the critical value of the elongation {x) = Xc + x, yields 
the Landau-De Gennes-like model 

Fc + Tl^tX + KoX^ - KlX^ + K2 x'^, (7) 

where the coefficients kq, ki and K 2 are considered pos¬ 
itive and is the tension of the RNA molecule at the 
critical value Xc of the elongation. The internal tension 
is therefore given by Tint — T^nt + 2 Koa; — + Ak, 2 X^. 

Upon substituting the last expression into Eq. an 
effective free-energy can in turn be defined that incorpo¬ 
rates the effect of the externally imposed tension r, that 
is 


F,~Fr 


- t) 




KiX 


^ 2 ^ 


( 8 ) 


The external tension tilts the internal free-energy land¬ 
scape Fx allowing transitions between the two possible 
configurations, folded and unfolded. The equation for 
the average elongation is now 


d{x) dFf. 

dt d{x) 


( 9 ) 


The effective free-energy is reported in experiments and 
simulations from which the values of the parameters and 
important biological information can be obtained [aisiia- 

mM- 

The numerical solutions of Eq. for the free-energies 
considered were obtained using the Euler-Maruyama 
scheme and presented in Figures 1 and 2. Critical os¬ 
cillations appear when the external and internal tensions 
balance each other out in such a way that the effective 
free energy F^ has two minima with similar energies. 
For sufficiently low or large external tensions, the free 
energy has a dominant minimum corresponding to the 
folded and unfolded states, respectively [see Figures la)- 
c)]. The noise induces critical oscillations between the 
two stable states separated by the free-energy barrier. 
The stochastic process associated with this oscillation of 
the molecule’s elongation is not a Gaussian white noise, 
as should occur for equilibrium fluctuations. On the 
contrary, it resembles a dichotomous stochastic process 
having a finite correlation time due to the intermittency 
which causes hopping events between both metastable 
states of the system. As it is well known, this inter¬ 
mittency, characterized by non-Gaussian distributions, 
is one of the characteristics of non-equilibrated processes 

EaUH]. 

It is also important to emphasize here that the critical 
oscillations associated to the hopping events cannot be 
understood as an example of an equilibrated chemical re¬ 
action 012]. The first aspect of the problem is to consider 
that unlike the concentrations of chemical species in equi¬ 
librium chemical reactions, a single molecule can never be 
in both states at the same time (coexistence). The second 



FIG. 2. (Color online) The entropy produced as a function of time 
(lines) by integrating Eq. for the single noise realization shown. 
Different intensities of the noise term cr = fcflTy where considered. 
Increasing noise intensity (larger cr) implies larger entropy produced 
because hopping events are promoted. Ip is the persistence length 
of the handles, L the length of the handles and x^ys the dimensions 
of the hairpin according to Ref. [2. 

aspect is that the thermodynamics of chemical reactions 
involves the concept of concentration, which in statisti¬ 
cal language entails the assumption of an ensemble of sys¬ 
tems. For dilute compounds, concentration and probabil¬ 
ity are related concepts. Thus, any chemical information 
and chemical kinetics description implicitly corresponds 
to an ensemble of systems, and not to individual ones. 
Our contention here, is that during a single experiment, 
each molecule cannot be in equilibrium because it jumps 
alternatively between the folded and unfolded states; its 
elongation is an oscillatory function of time. These hop¬ 
ping events imply a finite change of length of the molecule 
at a finite rate and thus cannot be seen as an adiabatic 
hopping process. It may happen, however, that in the 
case of a collectivity of RNA molecules the number of 
molecules passing from the folded to the unfolded states 
and viceversa, compensate each other, thus leading to the 
well known definition of the equilibrium constant in terms 
of the average reaction rates of a chemical reaction. The 
crucial point here is very well understood in statistical 
physics: Although the ensemble may be in a stationary 
state (there is not clear evidence of a thermodynamic 
equilibrium here), each system (molecule) develops a dy¬ 
namics, that is, a non-equilibrium process that should 
satisfy the well known fluctuation-dissipation theorem. 

The last result indicated in Figure Id) shows an ex¬ 
ample of the effect of a periodically time dependent force 
on the critical oscillations. Depending on the intensity 
of this oscillation, stochastic resonance and a more con¬ 
trolled jumps between folded and unfolded states may be 
obtained, thus allowing to infer more regular shapes for 
the free energy. It should be stressed that this oscillating 
force is a toy model of existent feedback mechanisms used 
in the experimental devices to maintain a constant ten¬ 
sion when spontaneous oscillations (folded and unfolded 
configurations) take place. Many details of this inter¬ 
esting aspect of experimental setups and their theoreti¬ 
cal modeling are well discussed in Refs. usiisoi for the 
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so-called passive or constant force modes. Here, we as¬ 
sumed a time dependent tension imposed on the system: 
t(<) = t[1 -I-dsin(tLit)], with d an amplitude factor and w 
the frequency, with the aim to emphasize that this kind 
of mechanisms may modify the response of the system 
as shown in Figure [^), leading to stochastic resonance 
events [S]. The results that can be obtained in this case 
have a direct influence on the probability distributions 
(or histograms) obtained from simulations and/or exper¬ 
iments, and therefore may influence some aspects of the 
determination of the free energy profiles. 



B. Entropy produced by a single molecule in a 
single realization 

For a single realization of the noise and, therefore, for a 
single sequence of hopping events of one RNA molecule, 
it is possible to calculate the entropy produced by nu¬ 
merically integrating Eq. Furthermore, a general 

analytic expression for the entropy produced can be ob¬ 
tained by realizing that the hopping events can be mod¬ 
eled with sum of Heaviside functions having an stochastic 
distribution of transition times U: 

x(i) = ^ AxQ[t — ti], (10) 


FIG. 3. (Color online) Three probability distributions as a func¬ 
tion of elongation for different times obtained numerically from Eq. 

a) Two probability distributions (modal non-Gaussian and 
bimodal) from the initial stage to the stationary regime in the case 
of a constant external tension, b) Two oscillatory behaviors of the 
bimodal distribution in the presence of a periodic external force, 
the alternative increase of probability in folded and unfolded states 
and the oscillation of the maximum corresponding to the unfolded 
state. 

of a probability density whose time evolution is concomi¬ 
tant with the Langevin-like equation Eq. (|^. 

For the present case with additive white noise, the 
Fokker-Planck equation for the probability distribution 
of X, P{x,t), associated to Eq. ([^ is 


with Aa; the elongation change between folded and un¬ 
folded states. The time derivative of this function is 
known as the comb function HI(t — ti) = ~ ^i) 

^^=Axm{t-u). ( 11 ) 

Using this result into and integrating over time one 
obtains that the entropy produced during a given elapsed 
time is 


'IT 

^^S=^Y.^T-TUt^)]■ ( 12 ) 

i 

For the realization of the noise shown in Figure 2, it is 
also shown the entropy produced in terms of time. The 
average behavior resembles a stair-like increase due to the 
hopping events between the two possible states. Thermal 
fluctuations in the unfolded and folded states of the dy¬ 
namics where subtracted in order to keep only the hop¬ 
ping events. 


III. ENSEMBLE DESCRIPTION 

The results obtained in the previous subsection refer to 
a single RNA molecule. However, the stochastic nature 
of the dynamics of the system allows one to introduce an 
interesting aspect of the problem that is associated with 
the collective behavior of an ensemble of equivalent sys¬ 
tems. This collective behavior can be analyzed in terms 


dP _ d 
dt dx 


dP 

innt -t)P + ksT— 
ox 


(13) 


For constant external tensions, it is well known the sta¬ 
tionary solution of this equation 

Peq{x) = (14) 

with Q = This expression pre¬ 

dicts a distribution function that may be bimodal de¬ 
pending on the balance between the external and inter¬ 
nal tensions, |26j . For example, for an externally ap¬ 
plied tension of 15pN [see Ref. 0], the value of the fit¬ 
ting parameter of the linear term of the free energy is: 

fc) 

An —T = —5.22pN. In this case, the effective free energy 
AFe has two well defined minima and it follows that the 
internal tension of the RNA molecule during the experi- 
ments was ~ 9.78pN. This is the tension associated 
to the handles at the critical elongation value, whereas 
the remaining 5.22pN are the apparent internal tension 
of the hairpin during the unfolding/folding process. 

In Figure 3 we show the spatio-temporal numerical so¬ 
lution of Eq. (131 for different shapes of the potential 


(also shown). For constant external tensions, stationary 
distributions are obtained that perfectly agree with Eq. 
(14). Their shape depends on the shape of the free en¬ 


ergy and can be bimodal or single-modal depending on 
the value of the external tension. 

However, in the case of stochastic resonance that in¬ 
volves the external time-dependent tension, T(t) = r[l -|- 
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dsin(wt)], the resulting probability densities are time de¬ 
pendent with no stationary limiting behavior. Oscilla¬ 
tions of each maxima of the distributions are clear and 
its amplitude depend on the amplitude of the externally 
applied force. For bimodal distributions there are two os¬ 
cillations coupled, that of the maxima with respect to an 
average value and that of the alternating change of the 
populations between the two possible states, folded and 
unfolded, see Figure 3. For single-modal distributions, 
only the first oscillation takes place. This is an impor¬ 
tant feature since if time averages are not taken over suf¬ 
ficiently large time windows, they can lead to a incorrect 
determination of a bimodal distribution and therefore to 
an incorrect determination of the free energy landscape. 

The probability distributions obtained in this subsec¬ 
tion show more readily the main features of the nature 
of the elongation’s fluctuations we already discussed in 
the previous subsection. For instance, if the external 
tension is above but near the bistable region of the force- 
extension curve, then the elongation of the RNA molecule 
is single valued (the free energy has only one minimum) 
but its distribution shows a large asymmetry that indi¬ 
cates a non-Gaussian distribution for the fluctuations of 
the molecule’s elongation. Such a situation is shown in 
the left panel of Figure 3 and is in good agreement with 
the histogram of Figure Ic) (not shown). Hence, at the 
edge of the bistable zone, these results clearly show a 
break of the classical condition for equilibrium fluctua¬ 
tions. 

This effect is more dramatic when considering the 
case when the tension takes values in the portion of the 
force-extension curve with negative slope. The distri¬ 
bution of these fluctuations is stationary and bimodal. 
The main implication of this result is that the collec¬ 
tive behavior reaches a stationary state characterized by 
a non-Gaussian distribution showing two peaks. These 
two peaks imply two most probable elongations of the 
molecule that do not correspond with the average elon¬ 
gation of the system, just at the center of the distribution. 
Hence, it is much less probable to find the system with 
the average elongation than in an extreme configuration. 
This fact contradicts the statistical condition of equilib¬ 
rium which is consistent with the condition of thermody¬ 
namic equilibrium: the average value of the elongation 
in the single peaked probability distribution is equal to 
the most probable value of the same variable m- 


A. Transition rates and mean first passage times 


Two important quantities characterizing the hopping 
dynamics of the RNA molecule between the folded an un¬ 
folded states are the transition rates and their inverse, the 
mean first passage times. In the present case, the tran¬ 
sition rates can be calculated using the Fokker-Planck 
equation (131 in the non-equilibrium stationary case |31j . 


As usual, the Fokker-Planck equation (13) is written in 
the form ^ = — Jyj, in which the current j is defined 


as 


J = 


ksT dx dx 


(15) 


where j is assumed to be a constant different from zero 
and D = ksT^. If we define Xf as the elongation of the 
folded state, then the complete solution of this equation 
in the interval [xf,x] is given by 

P{x) = (16) 


-F,{x)/ksT 

D 




Solving Eq. (16) for j and evaluating the result in x^ 
the elongation of the unfolded state, we obtain 

J ~ ^ _ p /i.„T . 7 . ~ jf 


r’”" gF^(z)/kBT 
Jxj 




(17) 

where jf and are the forward and backward cur¬ 
rents, respectively. For sufficiently large energy barriers, 
Fe,max » ksT, One may assume fast local equilibration 
of the system in each well and write P{xf) and P(x„) 
in terms of Eq. (14). In addition, the integral in the 
denominator of Eq. (17) can be calculated by extending 


the limits of the integration from [—oo, oo] and assum¬ 
ing that Pe(x) ~ Fe{xc) — ^F^^^{x — Xc)^, where Pg is 
the second derivative of Pe(x) evaluated at Xc, the lo¬ 
cal maximum of the free energy. These operations give 
the following expressions for the forward and backward 
currents 

1/2 


Jf = DQ 


f 


P. 


27vkBT 


^ — Fe{Xc)lkBT 


ju = DQ 


-1 


P„ 


-I 1/2 


2'KkBT 


^-F^{x^)/kBT 


(18) 


where Qf and Qu are the normalization factors of the 
folded and unfolded states, respectively. Einally, the 
transition rates are defined as kf = jf/nf and ku = 
ju/nu where 


nf — Qj 


= Q~n 


F. 


ej 


2'KkBT 


- 1/2 


,-F^{xf)lkBT 


P. 


2KkBT 


- 1/2 


,-F^(x^)/kBT 


(19) 


Here, we evaluated the normalization factors by assum¬ 
ing the usual quadratic approximation of the free en¬ 
ergy around each well, for instance: Pe(x) ~ Pe(x/) -I- 
iPg j(x — x/)^. Then, after substituting into the defini¬ 
tion of the rates we obtain the final expressions mm 


\jKj^e,c _r 

- e. L 


P’e(xc)-P’e(a3f) 


’‘1 = '’ ’ 2^ " 


ku = l 


^e.uFe c _ \ 

. [ fcflT J 


2k 


( 20 ) 
































7 


These quantities may be evaluated explicitly after fitting 
the free energy profiles reported in the literature [6] . Ta¬ 
ble 1 shows the values of the transition rates and the 


corresponding transition times evaluated for some cases. 
All these results fall within the same order of magnitude 
of the transition times reported in [^. 


Case 

X F [nm] 

Xu [nm] 

XTlnm] 

uW 

kf[l/s] 


ku[l/s] 

1 

1.93 

7.28 

3.71 

3.9 X 10-4 

2564.10 

5.9 X 10-4 

1694.91 

2 

1.84 

7.20 

3.85 

3.5 X 10-4 

2857.14 

4.4 X 10-4 

2272.72 

3 

1.96 

7.11 

3.92 

4.6 X 10-4 

2173.91 

4.5 X 10-4 

2222.2 

4 

2.13 

6.29 

3.6 

1.7 X 10-4 

5882.35 

11 X 10-4 

909.09 

5 

2.00 

7.41 

3.91 

3.9 X 10-4 

2564.10 

5.8 X 10-4 

1724.13 

6 

2.04 

6.90 

3.89 

4.9 X 10-4 

2040.81 

6.1 X 10-4 

1639.34 


Table I. Mean transition times and rates from the folded to unfolded states obtained from Eqs. 
energy profiles reported in [B]. 


(20) for the fitted 


CONCLUSIONS 

As a summary, in this work we have proposed an irre¬ 
versible thermodynamics of single-molecule experiments 
subject to external constraining forces of a mechanical 
nature. Our work is a powerful extension of Onsager’s 
formalism to the case of systems under non-equilibrium 
external constraints. Using this formalism we were able 
to calculate the entropy production of a small system in 
contact with a bath and to derive the general non-linear 
kinetic equations for the variables involved. A more fun¬ 
damental justification of our formalism can be given by 
postulating the existence of a partition function of the 
biological small-system in the presence of external con¬ 
straints and the validity of the thermodynamic stability 
criterion. Using these postulates, an analytic continu¬ 
ation of the free-energy into the unstable region allows 
one to relate the free-energy barrier with the entropy pro¬ 
duced by the single system during its time evolution. 

As an illustration, we have analyzed the folding¬ 
unfolding dynamics of RNA molecules subject to ten¬ 
sion. We showed that a single macromolecule evolves far 
from equilibrium executing critical oscillations when it is 
forced by the external constraints to remain in the un¬ 
stable region of its free energy. Since this free energy 
contains anharmonic contributions, these critical oscilla¬ 
tions are not a simple Gaussian process. This is a mani¬ 


festation of the unequilibrated nature of this process for 
which the entropy produced during hopping events was 
calculated. 

The approach we propose uses global thermodynamic 
variables for the single small-system in similar way as 
Onsager’s non-equilibrium thermodynamic theory. This 
fact enables us to derive non-linear differential kinetic 
equations governing the dynamics of the system in both, 
the deterministic and stochastic cases pn PS]. Unlike 
Onsager’s original theory, we have shown the validity of 
this description for systems far from equilibrium (non¬ 
linear regression in the presence of external forces) and 
not only near equilibrium without external forces (linear 
regression of spontaneous fluctuations). This opens many 
possibilities for the theoretical analysis and prediction of 
the behavior of single molecule experiments. 
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